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INTRODUCTION 


i 

For  most  of  the  tank  projectile  test  firings,  sky  screens  are  used  to  obtain 
the  flight  time  as  a  function  of  range.  These  range  data/can  be  used  to  extract 
muzzle  velocity,  velocity  decay,  and  drag  coef f icient .  '  The  most  common  technique 
to  extract  muzzle  velocity  is  to  use  the  difference  method  to  obtain  the  secant 

xi+l  -  xi 

velocities  (i.e.,  V—  =  — - The  first  few  secant  velocities  are  then 

x .  t  -  t 
1  l+l  i 

fitted  with  a  straight  line  or  a  parabola  to  obtain  the  muzzle  velocity.  There 
are  two  disadvantages  with  this  type  of  data  reduction  technique:  (l)  the  dif¬ 
ference  method  magnifies  the  error  in  obtaining  the  secant  velocities,  (2)  the 
velocity  profile  is  neither  a  straight  line  nor  a  parabola.  To  remove  these 
disadvantages,  a  closed  form  solution  to  the  equations  of  motion  is  used  to  fit 
the  flight  time  from  the  sky  screens.  A  closed  form  solution  is  possible  when 
the  drag  coefficient  is  assumed  to  be  a  linear  function  of  velocity.  Therefore, 
the  drag  coefficient  can  also  be  obtained  by  fitting  the  closed  form  solution  to 
the  sky  screen  data. 


ANALYSIS 


The  equation  of  motion  for  a  flat  fire  antitank  projectile  is 

MV  =  -*5pV“SCj)  (l) 

where 

M  ■  Projectile  mass 

V,  V  =  Velocity  and  time  rate  of  change  of  velocity 
p  =  Air  density 

S  =  Reference  area 

Cp  =  Drag  coefficient 

When  the  drag  coefficient  is  a  linear  function  of  velocity,  having  a  form  of 
Cq  =  Cp  (1  +  bV),  equation  l  becomes 

a 

M$  =  -lspVnCDa(l  +  bV) 
or 

V  =  -AV2(1  +  bV) 

(2) 

where  A  is  a  constant  and  is  defined  as  pSCpa,  and  b  is  a  slope  value  for  the 
linear  drag  expression.  In  most  cases,  b  has  a  negative  value. 


1 


The  closed  form  solution  for  velocity  as  a  function  of  range  (x)  can  be 

obtained  by  integrating  equation  2  with  V  =  dx.  The  solution  is 

at 


V  =  Vc (1  +  bV)e 


-Ax 


or 


V  = 


eAx  -  bVf 


(3) 


where  V 


I  +  bV,. 


with  the  initial  condition  of  V  =  V  at  x  =  0. 


Furthermore,  the  closed  form  solution  of  t(x)  can  be  obtained  from  differen¬ 
tial  equation  3.  It  has  the  solution  of 

t  =  tQ  +  -^-(e^  -  1)  -  bx  (4) 

with  the  initial  condition  of  t  =  t  at  x  =  0. 

Equation  4  expresses  the  projectile  flight  time  as  a  function  of  range  which 
is  in  the  form  to  fit  the  sky  screen  data. 

To  obtain  a  best  fit  in  the  sense  of  minimizing  the  squares  of  the  devia¬ 
tions,  a  residual  is  defined 


IT  ( t  • 

i=i  (  ic 


-  le ; 


2 


(5) 


where  t.  is  the  time  of  flight  calculated  from  equation  4  and  t^e  is  the  tyne  of 
flight  obtained  from  the  sky  screen  for  the  projectile  to  reach  the  ith  sky 
screen,  and  N  is  the  total  number  of  sky  screens  used  for  that  test. 


To  minimize  the  residual,  equation  5  is  differentiated  with  respect  to  the 
four  variables  (t0.  A,  VQ,  and  b)  individually,  and  then  set  to  zero.  They  are 


N 


N 

. (t ic 
i  =  l 


9 

i=i 


t  ie)  ^  ic  =  0 
3t() 

fje  =  0 

t  =  o 

3V0 


tie)^£  =  0 


(6) 


The  following  partial  derivatives  of  tjc  with  respect  to  the  four  variables 
can  be  obtained  by  differentiating  equation  4 


3ti 


?t0 


ic 


iA 

:>ti 


_L[1  +  eAxi (Axj  -  1)] 


AVC 


'X  \f 

o  vO 


i_(eAxi  -  1) 


av; 


ic 


(cAxi  -  1)  -  xt 


(7) 


With  such  a  complex  expression  of  the  four  variables  in  the  set  of  equation 
6,  the  easier  approach  in  solving  the  four  unknown  variables  is  by  using  the 
iterative  method.  To  obtain  the  corrections  for  the  variables  from  the  j  iter¬ 
ation  to  the  j  +  1st  iteration,  the  t*c  is  expanded  in  a  Taylor's  Series  and  the 
series  is  truncated  after  the  first  order  terms.  Therefore 


cic 


(j  +  1)  =  tlc(j)  +^i£|J(At0)  +lii£|j(AA)  +-^i£|j(AV0)  +  3tic| 


3t, 


(Ab) 


3A 


3b 


(8) 


s  t 

Equation  8  is  substituted  into  equation  6  for  calculating  the  j  +  1  cor¬ 
rections.  For  simplicity,  the  subscript  j  and  j  +  1  will  be  omitted  in  the  fol¬ 
lowing  expressions: 


*Rc 

at« 


:iRi 


iiV, 


ii 


-  t 


or 


+  (— -)?At0  + 
:)to  3to 


0tn  aA 


=  0 


N 

i=i 
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=  0  = 


.  I  (L  it:  L  icl 


atic  atir  at « c 

+  - —  — — 

lL  M 

ati 


aA 

at  i 


atr 


at 


1  £  x  7  „ 


at0  av. 


h/»tn  + 


at  ic  at 
a  a  av. 


^  ic„y 

+ 

3  tic 

»V0  '  0 

3b 

at  at  j  c 

^  1  /\\i 

3tlc.bi 

oa  avG  0 

3A 

ab 

+  (,'ti<  )-  ‘V  , 

+ 

3t  I,. 

'iVo 

)Vo 

ab 

w  o  o  o  o  o  c 

^  =  0  =  .5!,i(tlt.  -  tic>^ii£  +  2iici£iAt0  +  Ohs.  2ii£AA  +  ilii£  ilii£Av0  +  (i 

1  =  1  11  1C  3b  3t„  3b  3A  3b  3V(>  3b  ; 


i^)?Ab] 

3b 


(9) 


The 
form  of: 


set  of  equation  9  can  be  written  in  matrix  notation, 
M«  AB  =  C 


which  is  in  the 


3 


1 


where 


1  "W- 

3ti)2 

Z  3ti 

3ti 

S  3ti  3ti 

X  3tl  3ti 

1  »t0 

i=13t0 

3A 

1=1  3t0  3V0 

i=1  3t0  3b 

3ti  3ti 
l3t^ 

N  . 

n  ai 

i=l  3A 

i)2 

J  ati  3t± 
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N  , 
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3ti 
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3ti 

3V0 

N 

i=l 

n  atj.  »ti 

3ti  3t i 

N  3t i 

3ti 

N  3t t  3ti 

?  dii)2 
i-i  3b  j 

13t0  3b 

i=l3A 

3b 

illSVo  3b 

t> 

rr 

C 

N 

I  i- 1 (tie 

-  t  •  ) 3ti 

1L  at0 

AA 

N 
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i=l 

'  tic)ir 

AB 

AV0 
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N 

E  (tie 
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Ab  j 

N 

ii^l(tie 

(10) 


The  subscript  c  for  the  partial  derivatives  has  been  omitted.  There  is  no 
ambiguity  because  the  partial  derivatives  can  only  be  computed  from  the  closed 
form  solution.  The  corrections  are  obtained  by  inverting  matrix  M,  that  is 


AB  =  M_1C 


(ID 


The  iteration  will  be  terminated  when  a  prescribed  degree  of  convergence  or 
number  of  iterations  is  reached. 


DISCUSSION 


The  computer  program  listing  for  this  sky  screen  data  reduction  technique  is 
shown  in  appendix  A.  The  sky  screen  data  for  round  number  63  for  the  105-mm 
APFSDS  projectile  M8  33  fired  on  20  April  1984  at  Aberdeen  Proving  Ground  was 
chosen  as  a  test  case  for  this  data  reduction  technique.  The  computer  printout 
for  this  test  case  is  given  in  appendix  B.  The  results  show  that  it  only  took 
four  iterations  to  achieve  the  desired  convergence  even  though  the  initial 
guessed  muzzle  velocity  and  drag  coefficient  were  different  by  as  much  as  352 
from  their  best-fit  values.  Therefore,  for  reasonably  well-behaved  sky  screen 
data,  this  iterative  method  converges  very  rapidly. 


4 


% 


To  Investigate  the  deviation  between  the  path  the  projectile  travelled  and 
the  ground  distance  due  to  the  trajectory  curvature,  the  deviations  were  calcu¬ 
lated  for  ground  ranges  of  2,500  meters  and  3,000  meters  for  the  105-mm  HEAT 
projectile  M456  and  the  105-mm  APFSDS  projectile  M833.  The  super-elevation 
needed  to  reach  2,500  meters  for  the  M456  projectile  is  0.84  degree  and  to  reach 
3,000  meters  for  the  M833  projectile  is  0.41  degree.  This  deviation  is  approxi¬ 
mately  0.025  meters  for  the  M833  projectile  and  0.07  meters  for  the  M456  projec¬ 
tile  at  their  respective  ranges.  These  deviations  are  within  the  accuracy  of  the 
sky  screen  data,  therefore,  the  trajectory  curvature  can  be  neglected  in  the  data 
reduction. 


CONCLUSIONS 


The  only  limitation  of  this  sky  screen  data  reduction  technique  is  that  the 
drag  coefficient  has  to  be  a  linear  function  of  velocity.  From  results  of  some 
test  cases,  it  indicates  that  this  technique  still  rendered  excellent  results, 
even  with  the  drag  coefficient  deviating  slightly  from  the  linearity  requirement. 
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SYMBOLS 


A  Constant,  defined  as  -ji^'SCp^,  per  meter 

b  Slope  value  in  drag  coefficient  expression,  per  meter  per  second 

CD  Drag  coefficient 

^Da  Constant  in  drag  coefficient  expression 

M  Projectile  mass,  kg 

P  Air  density,  kg/meter3 

S  Projectile  reference  area,  meter2 

t  Time  of  flight,  seconds 

tc  Time  at  the  gun  muzzle,  seconds 

V  Projectile  velocity,  meters  per  second 

vo  Muzzle  velocity,  meters  per  second 

vc  Defined  as  - 2 - ,  meters  per  second 

1  +  bVQ 

x  Range,  meters 

(*)  Time  rate  of  change,  4— (  ) 

dt 

Subscript 

c  Values  computed  from  the  closed  form  solution 

e  Values  from  range  data 

i  Values  pertaining  to  the  i^  sky  screen 

j  Values  pertaining  to  the  jth  iteration 
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coimBft  mmm  listihg 


1 00 
1  10 
120 
130 
140 
150 
1 60 
1  70 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 


PROGRAM  MUZF  I  T  ( I NPU7  ,  TAPES-  INPU  1  ,  OU  T  PUT  ,  T  APG6-(j<J  I  PU  I  ' 
DIMENSION  X (16) ,TE<16> ,TC<16> ,TD(16) ,WXX(5) ,C<5,6> ,D(5) 

I)  1  MENS  ION  AA  (  4  )  ,  M  (  1 6 )  ,  X  X  <  1 6 )  ,  T  I  II  E  ( 8 )  ,  FF  1  4  > 

DIMENSION  R ( 71 ,0(7) 

DATA  XX/56.81 ,74. 81 ,92.81 ,110.8, 128. 81 ,239. 76,453, 22,779.9 3 
>982. 1,1251.18, 1479. O, 1741 . 45, 1984.35,2475. 85,2733. 21 ,2983.2 
DATA  AA / 3H I  0= , 3HA-  , 3HV0= , 3HB=  / 

DATA  F F / . 0 1 O , . 00005 , 30 . , . 00005/ 

DAT A  R/0 . O , 500 . O , 1 000. 0 , 1 500. 0 , 2000. 0 , 2500 . 0 , 3000. 0/ 

700  FORMAT  <8AI>  >> 

701  FORMAT <1611) 

702  FORMAT (8F 10.2) 

10  READ <5, 700)  TITLE 

IF (EOF <51  .NE.  0.0)  STOP 
READ ( 5 , 702)  (TE ( I ) ,1=1,16) 

READ  <5, 702)  (DID  ,  1  =  1  ,4) 

READ (5,701)  M 
L  -0 

DO  50  1=1,16 

I F"  ( M  (  I )  .  EQ.  0)  GO  TO  50 

L=L  +  1 

X  <L>  =XX ( I ) 

TE(L)=TE(I) 

50  CONTINUE 
NT =L. 

800  FORMAT <1H1> 

WRITE (6,800) 

600  FORMAT  < 1 HO / 1 OX , 8A 1 0 ) 

WRITE (6,600)  TITLE 

601  FORMAT ( / / 1 OX , *NUMBER  OF  SKY  SCREENS  =  *,I2) 

WRITE (6,601)  NT 

WR 1TE  <6, 80 1 ) 

801  FORMAT (//l OX, *RANGE  DATA-METERS* / ) 

WRITE (6,602)  ( X < I ) ,1=1, NT) 

602  FORMAT (  1 OX , 1 OF 1 O . 3 ) 

WRITE (6,802) 

802  F ORMAT ( / / 1 0 X , * T I ME  DATA-SECONDS*/ > 

WRITE (6,603)  (TE ( I > , 1=1 ,NT> 

605  FORMAT  <  1 OX , 1 OF 10.6) 

WRITE (6,803) 

803  FORMAT  (///10X,* INITIAL  GUESSES—*) 

WR I TE (6 , 604 )  (D<I),I=1,4) 

604  FORMAT (  /10X.F12.6  , F 1 2. 7 , F 1 2. 2 , F 1 2 . 6 ) 

NR=0 

E=0.0 

805  FORMAT  <//'4X  ,  13HITERAT  ION  NO.  ,  6X  ,  9HRLS I  DUALS  ,  7X  ,  1 4HPR0BABLL 
+// 1 OX , 43HC0RRECT I DNS  OF  CONSTANTS  FOR  EACH  ITERATION) 

WRITE (6,805) 


<• 


ERROR 


II 


b80  1  oo  SK-u .  o 

590  DO  130  J=l,5 

600  WX  X  ( J  )  o 

610  DU  130  1  1,4 

620  0 ( I , J ) =0 . O 

630  130  lUNTINUE 

640  TO- D ( 1 > 

650  A=D 1 2 > 

66C  VO" D ( 3 ) 

6  7 i.)  B~D  (  4  ) 

600  VC- VO/ ( 1 . 0+B*VO > 

690  A VO = A* VC 

700  AAVC=A*AVC 

710  VC2=VC**2 

720  002= VO* *2 


740 
750 
760 
770 
780 
790 
BOO 
8 1  O 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1 000 
1 0 1 0 
t.  020 
1 030 


DO  150  1=1, NT 
AX=A*X < 1 ) 

EAX=EXP ( AX  ) 

EAX 1 =EAX  -  1 . 0 

FTA= (-EAX 1+EAX*AX > /AAVC 

PTV=  -EAX 1  .  ( A*V02 ) 

PTB= (EAX1-AX) /A 

TC < I ) -T O+EAX 1 / AVC -B*X ( I  ) 

TD ( 1 ) =TE ( I >  -  TC ( X ) 

SR— SR+TD ( I ) **2 
C ( 1 , 1 ) =C ( 1 , 1 ) + 1 . O 
C (2, 1 )=C (2, 1 ) +PTA 
C  (3,1 )  =C  (3 , 1 )  i-PTV 
C ( 4 , 1 ) =C (4,1) +PTB 
C ( 1 , 2 ) =C ( 2 , 1 ) 

C ( 2 , 2 ) =C (2,2) +P  TA**2 
C ( 3 , 2 ) =C ( 3 , 2  * +PTA*PTV 
0(4,2) =0 ( 4 , 2 ) +PTA*PTB 
0(1, 3 ) =C (3,1) 

0(2,3) -C ( 3 , 2 ) 

C ( 3 , 3 ) =0 ( 3 , 3 ) +P  TV  **2 
0 (4,3) =0 ( 4 , 3) +PTV*PTB 
0 ( 1 , 4 ) =0 ( 4 , 1 ) 

C ( 2 , 4 ) =0 ( 4 , 2 ) 

0(3, 4) =0(4, 3) 

0 (4,4) =C (4,4) +FTB**2 
C ( 1 ,5)=C ( 1 ,5) +TD ( I ) 

0  <2, 5) =0(2, 5) +TD ( I ) *PTA 
C (3,5) =0 (3,5) +TD < I ) *PTV 
C (4 ,5) =C (4 , 5) +TD ( I ) *PTB 
150  OONTINUE 


1040  NR=NR+ 1 


1 050  E 1 =E 


<*■ 


t 


12 


f 


1 060  fc=0 . 6/45*SQRT ( SR/ (NT - 4 ) ) 

1070  WRITE (6,605)  NR , SR , E 

1 080  605  FORMA T ( / / 1 0 X , 1 2 , 2F20 . 8 ) 

1 090  IF  (ABS  (El --E )  - 1 .  OE -05 )  200 , 200 ,110 

1 1 00  1 1 O  1 F ( NR -15)  1 20 , 1 20 , 200 

1110  120  CALL  I NO (C ,  4 , 5  , WXX ) 

1120  DU  180  1=1,4 

1130  F=C(I,5) 

1140  I F ( ABS ( F )  .01.  F F  ( I > )  F=SIGN (FF ( 1 )  , F i 

1150  D  ( I > =D ( I ) +F 

1160  180  CONTINUE 

1170  WR I TE  < 6 , 606 >  (C < I , 5) , 1  =  1 , 4) 

1180  606  FOR’MAT  </10X,5E18. 6) 

1190  GO  TO  100 

1 200  200  CON  r I NUE 

1210  CALL  I NO (C , 4 , 4 , WXX ) 

1.220  DO  210  1  =  1,4 

1 230  210  W  X  X  < I ) =SQRT ( ABS ( C < I , I ) ) ) *E 

1240  806  FORMAT <//22X , 15HBEST-FIT  VALUES , 4X , 1 4HPRORABLE  ERROR ) 

1250  WRITE, (6,806) 

1260  DO  230  1=1,4 

1270  230  WRITE (6,607)  AA ( I ) , D ( I ) , WXX ( I ) 

1 280  607  FORMAT  < /  1 0  X , A3 , 5  X , 2E 1 8 . 6  > 

1290  804  FORMAT <//13X , 13HRANGE  -METERS  , 13HTIME - OBS  SEC  , 1 3HTIME- CAL  -SEC 

1300  *  1 2HT I ME-DIF-SEC/ / ) 

1310  608  FORMAT (1 OX, FI 3. 2, FI  3. 6, FI 3. 6, FI 3. 6) 

1 320  WR I TE ( 6 , 804 ) 

1330  DO  240  1=1, NT 

1340  WRITE (6,608)  X ( I )  , TE ( I )  , TC ( I )  , TD < I ) 

1350  240  CONTINUE 

1 360  V 1 =D ( 3  >  +D ( 1 ) *D ( 2  >  *D  <  3 ) **2  * ( 1 . 0+D ( 3 ) *D ( 4 ) ) 

1370  WRITE (6,609)  VI 

1380  609  FORMAT < // 10X , *MUZZLE  VELOCITY  AT  TIME  EQUAL  TO  ZERO  IS  *,F  10.2 

1390  V ( 1 ) =D (3) 

1 400  VC=D (3) / < 1 . 0+D (4) *D (3) ) 

1410  DO  250  1=2,7 

1420  VAX=VC*EXF' (-D (2) * ( 1-1 ) *500. 0) 

1430  V  < I ) =VAX / ( 1 . 0— D <  4 ) *VAX ) 

1440  250  CONTINUE 

1450  WRITE (6,610)  R 

1460  WRITE (6,611)  V 

1470  610  FORMAT </3X, 15HRANGE (METERS)  =,7F12.1) 

1480  611  FORMAT (5X , 15HVEL0CITY (MPS)  =,7F12. 1) 

1490  GO  TO  10 

1500  END 
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1510 

SUBROUTINE  INV (C , NC , NCS 1 ,WXX) 

1 520 

DIMENSION  C (5,6)  ,  WXX  <5)  , PIVOT (2)  .00 <5, 10) 

1 530 

NCT=NC*2 

1540 

NCP1=NC+1 

1  550 

DO  10  1=1, NC 

1 560 

DO  10  J  =  1 , NC 

1  570 

10 

CC(I,J)=C(I,J) 

1 580 

DO  20  1=1, NO 

1590 

DO  20  J  =NCP 1 , NOT 

1 600 

20 

CO ( I , J) =0. 0 

t  6 1  0 

DO  30  I = 1 , NC 

1620 

30 

CC(I ,NC+I)=1.0 

1  630 

DO  205  I = 1 , NC 

1640 

P I VOT < 1 > =CC < I , I > 

1  650 

DO  200  K=1,NC 

1 660 

PI VO  r (2) =CC  <K, I > 

1  6  'O 

1  26 

IF <K-I )  135, 130, 140 

1 680 

1  30 

DO  150  J=1 ,NCT 

1 690 

IF (PIVOT (1) )  134,210, 134 

1  700 

134 

CC (K , J ) =CC ( I , J ) /PI VOT ( 1 ) 

1710 

1 50 

CONTINUE 

1  720 

GO  ro  200 

1  730 

135 

DO  160  J =1 , NOT 

1740 

IF (PIVOT (1))  136,160,136 

l  75o 

1 36 

CC  < K , J ) =CC (K , J ) -CC ( I , J ) *PI VOT (2) /PIVOT ( 1 ) 

1  7  60 

1 60 

CONTINUE 

1770 

GO  TO  200 

1 7  80 

140 

DO  170  J=l,NCf 

170(1 

IF (PIVOT (2))  145,170,145 

1 800 

145 

CC ( K , J ) =CC ( K , J ) /PI VOT ( 2 ) -CC ( I , J ) 

1810 

1  70 

CONTINUE 

1820 

200 

CONTINUE 

1830 

205 

CONTINUE 

1 840 

GO  TO  250 

1 850 

2 1 0 

WRITE (6,600) 

1 860 

600 

FORMAT <//l OX, *DET  IS  EQUAL  ZERO*) 

1870 

250 

DO  300  1=1 ,NC 

1 880 

DO  300  J=1 , NC 

1  890 

300 

C  ( I , J) =CC(I , J+NC) 

1 900 

NCS=NCS1 

1910 

350 

T  F ( NCS-NC  >  500 , 500 , 400 

1 920 

400 

DO  420  1=1, NC 

1 930 

WXX  ( I  >  =C ( I , NCS) 

1 940 

420 

C(  I, NCS) =0.0 

1 950 

DO  450  1=1 ,NC 

1 960 

DO  450  J=1 ,NC 

1  970 

450 

C(I ,NCS)=C(I ,NOS)+C(I , J)*WXX  <J> 

1 980 

NCS=NCS- 1 

1 990 

GO  TO  350 

'2000 

500 

CONTINUE 

20 1 0 

RETURN 

2020 

END 
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